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ABSTRACT 


icmenecrs Soudies the roll loop of an operational 
surface to air missile as an example of a typical self 
adaptive control system. The missile studied also exhibits 
a phenomenum called ROLL WANDER,a low frequency oscillation 
in missile roll orientation. Because of this, the thesis is 
developed in three major areas of study. First, is the study 
Seeeone missile roll control loop. Second, is a study of the 
phenomenon of ROLL WANDER,with the production of this pheno- 
menon in the computer simulation. Third and finally is the 
implementation, in the computer simulation, of a self adaptive 


notch filter to suppress the occurence of ROLL WANDER. 
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ieeecOMUROL SYSTEMS 


A. INTRODUCTION 

The objective of this thesis was to study a self- 
adaptive control system through the use of a computer 
memileatciton. ithe control system studied was the roll control 
loop of an operational surface to air missile. 

Hoes Darutcular control loop uses a bistable actuator 
(a nonlinear element) as the primary mover of the loop. 
This method of control, a full throw "bang bang" system is 
wommonly used in control systems. The control loop is driven 
either full positive or full negative, depending upon the 
polarity of the error signal. Since the system reacts with 
full response, the output is driven rapidly in opposition to 
wm@e error and the error signal is eventually driven to the 
Opposite polarity. The control then changes sign and the 
output is driven in the opposite direction. In this manner, 
eeemerror signal is nulled out. 

Bang bang control systems are in wide-spread use for such 
Py cbens as surface to surface missiles, surface to air 
missiles, and some types of surface and submarine torpedoes. 
All of these systems require an accurate, reliable type one 
system producable at a minimum of expense. A full throw 
control system can provide these requirements. 

imeune case Of a surface to air or air to air missile, 
however, there are added problems that the control system 


must be able to cope with. It must be able to operate at 





Vervous aleicudes and at various speeds. The correct control 
response to an error signal of a given magnitude at high 
altitude (thin atmosphere) would be a drastic over response 
to the same error signal at low altitude (dense atmosphere). 
memeaorrecry 1O0r this, the control system must be self adaptive. 
meemuss respond to the error signal at a rate appropriate 
with its medium. Another constraint placed upon the control 
system is the requirement not to disturb the missile guidance 
and seeking/tracking mechanisms. While a full throw control 
system is inexpensive and reliable, it does introduce a 
Semsuaie Oscillation in the position output. This can be 
disruptive to the other missile systems mentioned. 
Thessurtace Go air guided missile studied in this thesis 
Delt zes a bane bang type of control system which is both 


self-adaptive and non-disruptive to other missile functions. 


Eee OLB RING IN THE ROLL SYSTEM 

pPigtweoml—| 15 a block diagprameof the roll control loop. 
The bistable element provides the control command to the 
PomMmMcemeroels surtace actuators. The use of this nonlinear 
mlenent in the control loop causes a iimit cycle oscillation 
Domo pcre btn Lom pnts Jamit cycle OSclllation which provides 
the self adaptive feature of the autopilot, and at the same 
time allows for smooth loop response while using a bang bang 
tvwoeecomenrol system. With no error or sorand signal present, 
the limit cycle oscillation provides a constant amplitude 


dither signal to the tail control surface servos. The dither 
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Wim Zero error signal has equal negative and positive 
Semeoeonenvs;, and as a result, the clockwise and counter 
el@ekWwise roll torques produced by the tail surfaces cancel. 
faewmrssile therevore dithers slightly in roll but maintains 
MEcOlGcale average roll orientation. The amplitude of the 
dither signal, and thus the adaptive gain characteristics of 
tee auvOpilot 18 a@efiunction of G, the roll moment coeffi- 
cient. G is a measure of the aerodynamic gain, a property 


of the medium in which the missile is traveling 


C. SELF-ADAPTIVE FEATURE 

The autopilot derives its self-adaptive feature by 
Tears tne amplitude of limit cycle oscillation, or dither, 
in the roll control loop. A study of Figure 1-1 shows that 
There are three basic loops. The outer two loops provide 
position (¢') and velocity (¢') feedback paths for stability 
and damping. The inner loop also provides some rate (6!) 
feedback, but of greater importance, the inner loop coupled 
with the bistable. element produces the limit cycle oscillation. 
ic roll Loop of the missile was used for the dithering loop, 
Smcemeanhieriae if the other control loops would not be 
meaccucel, Dithering in the pitch Loop would be intolerable 
due to the amount of body motion sensed by the head rate 
gyros. Similarly, a dither in the squeeze mode would not 
provide the aerodynamic adaptive gain features. 

Hoawet@om 1 is the dynamic equation for the rigid body 


motion of the missile. 


lk 





b: G3+FO ) 


ae 0 eG: (la) 
2 (s-F) 


where: ( the aerodynamic moment 


{| = the aerodynamic damping 


This transfer function, called the "aerodynamic roll . 
transfer function" produces the inherent self-adaptive 


feature. 


G= 1481\Sd M'57.3 C,. 


dl 
= (1b) 
where -\ = static air pressure ratio 
) = mach no. 
Te =) Fo iene reade 
Cis = POimecanemae. COet Ti erent 
2. c 
f= 43 | A\Sd my Seo Cir) (ale) 
a | oe 
where: \/ = velocity 


5 db 


aerodynamic constants 


As can be seen by Equations (1b) and (lc), the gain of 
the roll transfer function is dependent on flight conditions. 


Therefore, as the missile's environment changes, so does the 


an 





pommeol the roll transfer function and, accordingly, the 
mipeeewuide OF the position loop's limit cycle oscillation, or 


Gasgher . 


Pe ROGL CONTROL 

The controlling action of the self-adaptive loop is 
shown in Figure l-2. The tail control surface drive motors 
eee aoriven by the output of the bistable element at a constant 
fiawe in C€ither the positive or negative direction. When an 
Peror Signal 1S present, it appears as a de bias signal 
ened With the roll loop's natural dither oscillation 
(Figure l-2a). The algebraic sum of the two signals is fed 
to the bistable element (Figure 1-2b). As its name implies, 
the bistable element has only two states, full positive and 
full negative. When the value of KF. (the summation signal) 
passes through zero the bistable element shifts states. The 
fom mom Ol the dc bias produced by the error signal is to 
cause the square wave output of the bistable element to 
remain positive or negative longer during a complete cycle, 
as appropriate (Figure (l-2c). Since the tail control 
SuUriace motors act_to integrate this output, the result is 
a dithering movement of the control surfaces modulated by 
the appropriate control surface movement to null the roll 
errors tenal. | 

As previously discussed, the gain of the roll transfer 
function determines the amplitude of the dither oscillation. 


As a missile flies higher, the gain decreases and therefore 
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the amplitude of the dither decreases. The effect of the 
dec bias is more pronounced with the result that the control 
loop is more sensitive and has higher gain. Conversely, 

at lower altitudes or higher speeds, the dither amplitude 
memereaver and the error signal has less effect. The loop 
gain is reduced (Figure 1-3a). 

The dither imposes a maximum signal strength allowable 
for the error signal. If at a given dither amplitude an 
error signal is encountered which is greater than the dither 
emplitude, then the oscillation is SeEeea completely to 
one side of the axis (Figure 1-3b). The bistable element 
stays continuously in one state and the tail control surfaces 
are driven continually in one direction. Loop response 
diverges and the loop becomes unstable (Figure 1-3c). To 
prevent this occurrence, the error signal is electronically 
limited. 

pieeculm response in the roll loop is obtained by super- 
Ompesing the de error signal on the dither. In fact the 
error signal is Neto pllerdcestenal. buvsais also a time 
varying waveform. It is important to remember, therefore, 
that the dither frequency must be sufficiently higher than 
the missile control frequencies (error signal frequencies) 
so that the gain adaptive feature can be maintained. This 
imposes a severe lower bound on an acceptable dither 


frequency. 
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E. DESCRIPTION OF LOOP COMPONENTS 

Figure 1-1 can be divided into ehree Major pOrvuions : 
exes tall control surface actuators, the air frame, and the 
autopilot electronics package with associated sensors. 

The original design work on the missile started with a 
Biven set of actuators and the air frame's equations of 
emronw tne autopilot control equations were then written 
to drive the air frame with the given hardware. 

Three major parameters are input or sensed within the 
roll convrol system. They are oo? the commanded roll angle, 
K Fo, the bistable element forcing function (commonly called 
"delta dot called for") and $', the tail roll rate as 
measured by the roll rate gyro. o! TSemvecraved Dy tne roll 
free gyro to produce the sensed missile roll angle ¢'. . 
and o' are then differenced to obtain o> Che role error 
Signal. 

Pouce romuOrcsne Hn uUneCl aon 

Equation 2 is the equation of the basic roll forcing 


minGu1on lnbouwm@e bistable switcher. 


/ | a \ 1 / 
Xs Fs = (g- a.) Lim + - ye \K<é < AsiSC+ 85) O | (2) 
By i ee) 


Wie iee! K« Fs =e ecemunrad Signal 


the electronically limited error signal. 
If the error signal is allowed to get 
too large, the loop will stop dithering 
and diverge 


& 
= 
r 


18 





/ 
0 =erohiueirace OUGDUL OL the roll rate gyro 


Ko, Xo are system constants 


Preliminary analysis showed that high frequency 
resonances caused by the missile's elastic body bending mode 
were possible. To suppress signals at resonant frequencies 
and at the same time maintain the dither frequency, a 


cascaded set of high frequency (Wpsiter > Was ther? notch 
filters were added. The final form of the roll forcing 


function (Equation 2) becomes: 


| oe ~7>.\A! 
Ks ts = GU Oeban + 231( Ker 5188) Ga) (3) 


(1+ sa65 





where 


Ge = = — +/ 5S \2” als 
| sas (0) | 
CG, = Lis LOD « (aS. )* 


es ai 
ji+ 460 S+(z55) | 


Gp = [ls BaF 9 +6 ow 5: ee 





an DLs (95 a ‘ 
I+ Gab 3 +(2B5) [| + Bes + S2V] 
[1 + 2608) s + (Be) *] - 
ro 
[1+ 26s tee)" J 


Ie, 





eee Rate Gyro 


ie wbmaMeater fuUNnectiOnesrom the roll rate gyro is 


a 
| 
Q' ieee 2) 
-/ 
where: 0) Seine) GOmoumOuG ot ema 1 


missile roll rate 


& 
I 


mechanical damping coefficient 


Uw) 
il 


Be Ole Pree Gyro 
Thesvoblwiree ayro has a transfer MUNG tOre Ons lita UY 
and 1S considered linear over all normal operating ranges. 
4, Tail Servo Transfer Function 
inemucilaeServo  transter Lunction (using a rigid body 
model vice an elastic body model) can be described as 


x GES 
ODS S 
x. e G5) 


where: fie GhevactWator timemdelay 
=e approximation to the extra phase lag intro- 
duced by high frequency noise filtering, switch 


hystersis and aerodynamic plane lag due to the 
Vor oLtherimae 


Pee oUetace Lance. Moment 


Shown as a positive feedback path around the tail 
servo transfer function is K Kn: Ro RepEesenvs che lance 
moment due to tail position. Aerodynamic loading forces 
alevulmewon tne tail control surfaces produce this hinge move- 


(temomivenageuene OQULPUG Of the tall actuators. The loop 
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| 
: 


is shown as a positive feedback loop, since the greater the 
Sauective tail angle, the greater the hinge moment torque 
produced; which, in turn, tends to increase the effective 
tail angle. Kan is a conversion factor, and is given as a 
Svyecvem constant. 

In this thesis, the effect of hinge moment relative 
EorOuner System parameters was negligible; and in all 
Simulations it was set to zero. Hinge moment was included 


mie une model for completeness. 


6. F_ and F 
a squeeze 


Fs mono up Len COMuro! sipnal. It is produced in 
Pac DiGCh comerol Loop, then summed with the roll control 
signal as an input into the bistable element. F squeeze 
the control signal generated by the autopilot to compensate 
for the misalignment of the tail surfaces. Misalignment, 
resulting from mechanical tolerances and system slop in each 
Pec outm vail SUrtaces ‘Causes Lorques on the missile 
frame; counteracting Gorguesuileav wou proauce roll or pitch, 
Mid improduee frame Stress or squeeze. 

7. Biases 

The effect of biases in the roll loop is to cause 
Poemmussile to assume a roll attitude other than that estab— 
lished at time of launch. Maintenance of roll orientation 
iemeecential to the guidance portion of the missile. The 
various components of the roll loops in addition to the pitch 
GOOD en rougn Fe and the squeeze Ssignai, F all 


squeeze’ 


introduce biases. 


aul 





miese Loop biases are counteracted by the filter 


peock: 
E(s) = a. (6) 


The equivalent roll error of each contributor can be deter- 
fegeceDy SOlving for 6 in the following simplified static 
roll control equation (7), with each bias substituted in 


terms of deg or deg/sec in the appropriate term. 


(7) 


B= k U8) 0, 33iko 4 Keh)] + Sec 
S Ko 


eZ 





i. ewe CI NG Missa Le 


A. INTRODUCTION 

There were two major areas of study in the research. 
fmist was the construction and testing of an accurate 
computer model to simulate the missile. Second was an 
investigation of an abnormal phenomenon experienced in the 
missile roll control loop called "roll wander." An explana- 
tion of this problem and a proposed solution will be 


Meeoented in subsequent chapters. 


B. SIMULATION 

Due to the complexities of the missile control system 
Mmoemrdeoetnercent nonlinear operation introduced by the 
Distable element, it was decided to model the control system 
using a completely digital simulation. ie model was there- 
fore constructed using IBM's Digital Simulation Language 
(DSL). DSL was chosen for the model because under its format 
complex transfer functions could be described by single 
blocks with accompanying data statements. This eliminated 
Tne need for writing state variable equations for as high as 
eighth order systems, and thereby significantly reduced the 
amount of coding required to model the missile. A listing 
of the DSL program with a more detailed commentary on the 
operation of the DSL program is provided in Appendices A 


and B. 


= 





Pee ene neeure 2-l is the missile roll control loop as 
it was arranged for encoding into DSL. The letter names 
for the different function blocks correspond to the variable 
memes Used in the simulation progr:m. The DSL function TRNIR 
used to represent the various transfer functions imposes one 
Mmemerrecton. tt requires that the polynomial in the denomi- 
nator be of a greater order than the polynomial in the 
mera bor. the reason for this is because the TRNFR block 
determines the correct output response by solving the 
Pomreosoponding system of first order differential equations. 
This requires the polynomial in the denominator to have the 
meee s. order to allow for the construction of the correct 
ferst order system. To accommodate this restriction, certain 
blocks (B and HP) were formed by combining other blocks (as 
seen by a comparison of Figure 1-1 and Figure 2-1). This 
altered the appearance of the block diagram, but not the 
Pemev1Or of the loop. 

Block D represents the mechanical system's response time 
to the bistable element's output signal. A finite time 
occurs between the electrical command and the mechanical 
movement of the tail control surfaces. This was modeled in 
the computer program by the DSL function block DELAY. DELAY 
is a dead-time function. Sample points of the input expression 
ereemoavea at specified intervals. The output is set equal to 
zero until the delay interval expires at which time the output 
has the value of the input at the start of the delay. If 


necessary, linear interpolation is used between sample points. 
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tes aoility, gained using an all digital simulation to 
meeurately and easily adjust the delay interval, was one of 
weenie Or Tactors in choosing the digital simulation over an 
equivalent analog simulation. 

The bistable element was modeled using a standard 
PORTRAN logical IF statement. The bistable element has no 
gor polnu, Only Tull positive and full negative. The 
logical IF was written so that when K,F, was greater than 


6° 6 
zero, the bistable output was positive. When K 


sls Was Zero 
or negative, the bistable output was negative. 

The aerodynamic gain constants G and F were calculated 
at the beginning of each run. The parameters necessary to 
calculate these values and the parameters specifying the 
Other systems constants were specified at the beginning of 
peerver run. These constants could be changed during or at 
eaewend Of the run, thus giving the capability of multiple 
mois within one DSL job. 

Dirttieulty was initially encountered due to the step size 
Peoed. = Ine TRNFR blocks are comprised of a series of inte- 
erators simulated by DSL by a call to an integration routine 
of a Psp celrica by the user. To conserve computer time, 
meEeenca inverval fourth order Runga Kutta method was used. 
However, it is a4 property of this integration routine as 
well as many other numerical integration routines, that if 
the frequency components of the function to be integrated 
approach the frequency of the integration routine, specified 


Dyeeme swceGUtLine Step size, the integration method produces a 


26 


a 





divergent error and becomes unstable. Such was the case with 
Petomocoegram. initially, a step size of 0.005 was used. 
meomhems were encountered with the inclusion of the high 
/meeemency notch filters. The notch filter represented by a 


Genominator of: 


Ze 
moo.) Ss + (_s (8) 
1000 {Q00 
has a center frequency of 159 hz. 
oMebarivy, the noteh filter represented by: 


630 


1+ 20(0.9 ae 
29005 + (8) il 


has a center frequency of 108 hz. Using a step size with 
feces ration frequency of only 200 hz was insufficient, 
em@e@wcaused divergence. 

Optimal step size was obtained by testing the loop using 
a variable step size Runga Kutta method with different error 
limits specified. Dy spectnyvame agzolerance of 0,001 
absolute error an average integration step size of 0.0004 
was obtained. UWpisnmeamurcateney Of 2500 hz was sufficiently 
high to avoid numerical integration problems. During the 


actual runs, a final step size of 0.0001 was used. 


fe eoliNG THE MODEL 
Swope Or vine problemot determining the program inte- 


Preeronesctep Size, the only difficulties encountered with 


Zi 





Maemanacrtal Simulation runs were those of writing the standard 
meOocram LOfzic and debugging routine syntactical errors. Once 
memouvructed, the model was tested under normal operating 
@emareltons to ensure Chat it behaved the same as the real 
missile. The model was allowed to "fly" with no input 
(Figure 2-2). The initial jump noted is the transient 
behavior of the loop. Since the parameters of the loop were 
fenvinualiy being changed from one run to the next, initial 
@emarcLtons Of the integrators were not calculated; rather all 
mmeeeravors were set initially to zero, and the transient 
menmweonce accepted. ‘Ihe low amplitude oscillation is the 
bistable dithering action. The dither appeared in the 
mesowvren loop because of the step size of integration used. 
Eye aecreasing step size, the amplitude of the dither could 
Bemreduced. However, decreasing step size proportionately 
increased run time. A compromise was reached using a step 
paeze Of 0.0001. 

iaemerased 1oop) transter function for a roll command 
(¢/$.) peonned by writing the loop equation for Figure 1-1, 


(with K = 0) is: 


15, ® 
O° Kgs? 1+ 25 +fa)* Kokske (88 ) 


10 


(10) 


331G8| ker SEG, ea kke (BI ) 


<000 
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where: Ge 


G65 =SsrOrsrenal motch filter 


tHareorenT tune tton of the roll rate gyro 


eee ene vne iinal value theorem for a step input 


(¢, = $./s), to Equation (10) gives: 


io de & Ke Ke Ce) 
ie: eke {2\(16 3 ) Gab 


Similarly, the closed loop expression for a tail induced 


torque (G75 ,)5 also obtained from Figure 1-1 is: 
p= 8d 
Se Se ene 
2(.3) oom 
7 -ky (2 ie eats ) | aoe 


Kos (t+ 262s +(an3) ) + Ko Kske (2) 


— ae eee 


2-000 


cre 


ey 





iene 2Or a step input (6, pcos, Ss)... and applying the final 


weime theorem to Equation (12) 


eS on 


The model was thus tested using a step input of -45° for 
>, and Pescep input of 3° for bo. The response to the step 
mapuc of b. is shown in Figure (2-3). For a roll command or 
mero! error signal, the missile will roll and maintain the 
desired angle. The response of the model compared favorably 
with that of an actual missile. ‘The measured time constant 
of the missile was approximately 0O.21 seconds over the various 
flight conditions. The measured response of the model was 
wero seconds over the same flight conditions. 

The transient response time to a step torque, equivalent 
Mora step input (6 .) of 3° was approximately 0.5 seconds for 
the actual missile and 0.62 seconds for the model (shown in 
Figure 2-5). 

fests were soneeesed on the missile model under two sets 
oni flight Congres Ons «covering the limits of normal operation. 
Palen & Ponaiticn No. 15 represented normal altitude and 
velocity, whereas flight condition No. 24 represented low 
alvacude and normal velocity. The self-adaptive nature of 
the control loop can be seen in Figures 2-4 through 2-10. 

The amplitude of Keke changes according to the aerodynamic 


conditions, but the response (o/, or o/6,) remains constant. 
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Prous 2-2 


¢' vs TIME with No Inputs 
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BIGURE 2-3 


¢' vs TIME with ~45° step input @ t=3.0 sec 
Geeticwe conaition Fi5) 
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FIGURE 2-4 


K F vs TIME with -45° step input @ t=3.0 sec. 
(°¢fight condition #15 ) 
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FIGURE 2-5 


$4 vs TIME with +3° step torque (6) at t=3.0 
@elicht. Condition 15 ) 
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MiEGURE 2-6 


K 5¥ VeeBiME with +3° step torque (5. ) at t=3.0 sec, 
( Sefight Conareron #15 >) 
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BEGURE, 2-7 


64 vs TIME with a -45° step input at t=3.0 sec. 
Gerlvgit condition #24 ) 
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BeGouns 2-8 


Zoe iMeE wrth a —45° step input at t=3.0 sec. 
i Pian condition #24 ) 
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PREURE 2-9 


¢'vs TIME with a +3° step torque (6) 
(emieeght Condition #24 ) 
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PacuRe 2-10 
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(°¢fight condition #24 ) 
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D. CONCLUSIONS 

The simulation flights shown in Figure 2-e through 
Figure 2-10 were compared to the results of test flights 
Periormed by the civilian contractor. On the basis of these 
comparisons, the conclusion was reached that the model was an 
Meerrave representation of the actual missile control system 
and that the model could therefore be used in the study of 


the roll wander phenomenon. 
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Tift. ROLL WANDER 


A. INTRODUCTION 

As stated earlier, the two major areas of research were 
the construction and testing of an accurate computer simula- 
tion for the missile and the study of the phenomenon of "roll 
wander." 

Sais esundesired rol] loop excitation or roll wander re- 
wees 21rom noise originating in the tail control surface 
@epuabvors. This noise is sensed by the roll rate gyro and 
is thus introduced by the roll rate gyro into the rate feed-— 
Beek loop as a roll rate error signal. The production of 


roll wander then occurs in one of two fashions. 


B. LOW AMPLITUDE ROLL WANDER 

The noise signal introduced is enhanced by the effects 
of the roll loop nonlinearities. The effect of these non- 
itjmearities is to generate harmonics of the dither frequency. 
ne Se harmonies in turn BOG Waki eu EUndamental- Gither: 1re— 
quency to produce a distorted wave form of a frequency equal 
to the Poneence! frequency but with zero crossings within 
one cycle unevenly spaced (Figure 3-1). Since the bistable 
element is triggered only by these zero crossings, the average 
frequency of the bistable remains the same; but the time in 
its two states becomes unbalanced (Figure 3-2). The neutral 
position of the tail surfaces is the integration of the 


bistable element output. Therefore, the tail control surfaces 
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PAGuURE 3-1 


VouliMe Showing Harmonic Addition 
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PEGuURE 3—Z 


Bistable Output 
(Seearmonics Commence at t=1.0 sec. } 
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begin to move off of their zero positions and the missile 
feels oance no roll command was introduced into the loop, 
Seeesultant dc error signal is generated in the roll position 
loop (outer loop). This dc signal causes the input to the 
bistable element to be shifted accordingly above or below 
MiiemeaetO axis. oince the distortion of the waveform is 
@avoead Dy harmonic addition, and the fundamental frequency 
Seeene bistable remains unchanged, a shifted position exists 
for the waveform where the zero crossings once more become 
ieee tne result is that the missile initially rolls off 
Pee on antroduction of the harmonic, then rolls back 
(Figure 3-3). If the phase of the harmonic wave is shifted 
rapidly enough to repeatedly offset the resultant waveforn, 
Mmgemmssile will begin to wander in roll at a frequency of 
1.2 to 2.5 hz trying to compensate (Figure 3-4). In the 
faouaile missile, uncorrelated white noise from the tail con- 
wero! actuators produces this phenomenon. The broad band 
noise contains components at the harmonic frequencies of the 
@ecaer Signal. Since this noise is random uncorrelated 
noise, inere Exists no coherence in the phase of the various 
spectral PonoCnente. Therefore, the action of the noise 

is the same as the addition of the fundamental dither fre- 
quency with a series of harmonic frequencies shifting 
rapidly in phase. Noise with very low amplitude, only 10% 
of the amplitude of the signal (SN = +20db) is sufficient 


to cause the onset of roll wander (Figure 3-5). 
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mEGURE 3-3 


vce uilibecwovaing themptfect of a Fixed Phase Harmonic 
( Harmonic Commences at t=1.0 sec. ) 
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FIGURE 3-4 


¢'vs TIME Showing the Effect of a Varying Phase Harmonic 
(Harmonie Commences at t=1.0 sec. ) 
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cu ies—-5 


¢' vs TIME Showing the Effect of Adding White Noise 
(ecenstcant level ) Noise Commences at t=0./ sec. 
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Note: Initial excursion is due to the initial conditions 
used in the simulation. 
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C. LARGE AMPLITUDE ROLL WANDER 

As the noise amplitude increases, the amplitude of the 
wander increases, since the harmonic components of the noise 
femal also increase. The frequency of the wander, however, 
remains constant (Figure 3-6). The action of the noise is 
to reduce the low frequency roll loop gain. The decrease 
in the low frequency roll loop gain will decrease the roll 
mereemedamping in the closed loop response, resulting in 
maeser wander excursions. If the noise gets large enough, 
the bistable switch gain may be decreased to where the roll 
system will enter into a sustained low frequency oscillation 
of extremely large magnitude. The dithering action of the 
loop is stopped and the entire system operates in a saturated 


limit cycle mode of operation (Figure 3-7). 


D. SIMULATING THE ROLL WANDER 

1. Noise Source 

One of the more difficult problems was that of accu- 

rately aiimtlasiine wae roll wander phenomenon. Analog simula- 
tions performed bythe civilian contractor had been performed 
using 400 ae Epand limited gaussian noise. Toproduce this 
digitally, a FORTRAN subroutine was written utilizing a 
random number generator and a digitally simulated 400 hz low 
pass filter to generate the appropriate band limited noise. 
mic daigital simulation of the transfer function of a low 
pass filter requires the integration of all data points to 


Peoguce the filtered output. Therefore, to preclude the 
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PIGURE 3=6 


mvs Time Showing the Effect of Adding White Noise 
meinearly wncreasing level ) Noise Commences at t=0.7 sec. 
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REGU RE 3 —7/ 


onl ME onovonc LogegIinstablity due to Saturating 
the Loop with High level White Noise 
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fMecesosity Of those time consuming integrations during the 
actual missile simulation, twenty thousand points of 400 hz 
memalimated white noise were computed, then written on disk. 
A subroutine in the DSL program called NOISE was called each 
meme increment. With each call to subroutine NOISE a read 
was performed from disk and the required white noise was 
muepiaed as output. Since a DELT (integration step size) of 
0.0001 was used in the simulation, two seconds of white 
noise were available. At the end of that time, the noise 
record was rewound and reused. The effect of this repetition 
meine NOolsSe record every two seconds of simulation time on 
Memes randomness of the white noise was considered negligible. 
eam le Of the noise used in producing: the wander shown in 
Figure 3-6 is shown in Figure 3-8. 
fwaeeuarmonie Generavor 

Since low amplitude wander vice large amplitude roll 
wander was of prime importance, it was decided, that to 
effectively study the phenomenon with the hope of preventing 
it, a pure poe nerte waveform rather than broad band noise 
should be used as an I gyore ge 

An eonacer second harmonic could be generated by 


using a square law device. 


then if: BI = BI Cos. 3 
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FIGURE 3-8 . 
White Noise (linearly increasing level ) 
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HARM = BI COS” 9 
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Bf = BI + HARM 


isl == jee COs ot Bie 
ma 


4 3 Cog26 
ae (3 + 3 COS2eG) (15) 


x 
Mmremiarmonic term however carries with it a de term which 
memeed cause a roll offset, and more importantly, might hide 
miemac Crror term generated by the control loop. To eliminate 
mse OC level, the output HARM1 was passed through a three 
stage low pass filter with a -3db frequency of 10 hz and 

eweto db/oct roll off. 

As discussed earlier in the chapter, when a roll 
meet CXisSts, adc error term is produced to null out the 
loop. When the harmonic is added to the loop, an immediate 
de term is produced moving the fundamental waveform off the 
mee axis. Since the harmonic is obtained by squaring the 
ememat in the outer lec. firs signal has such a de error 


term then: 


HARM = BIS 
3 tse. 4 om ae 
O ii 
= BIT + 2BI *BL COSsa 6 
+ BI, (s+5 EGsc6) 
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Equation (16) represents a distorted waveform. As 
BI. (the de component) increases, the effect of the harmonic 
memponenvt decreases. To avoid this problem and again to 
enable viewing only the wander phenomenon, the waveform used 
momproduce the harmonic was passed through another three 
moose low pass filter similar to the one described above, 
(Figure 3-9). Thus at anytime, only a pure sinsoid at a 
harmonic frequency was introduced. 

CUMUtOMmwOmnbitcepmase locked harmonic is shown in 
feeeure 4-2. The excursion caused by the addition of a har- 
memie Signal producing a signal to noise (harmonic) ratio 
of approximately +12 db was -11.18°. To produce sustained 
roll wander, the phase of the harmonic signal had to be 
Pommucad rapt aly 4n relation to the fundamental. To accom- 
Dplish this, a DSL DELAY function block was used. The amount 


ef the delay was varied by setting it equal to: 


Shift = ¢ Sin 24.1 
: ve) 
where: 
+, = 10 hz 
ty = fundamental dither frequency 
then: HARM'(s) = HARM (S) e7(shift)s 


The resulting waveform is shown in Figure 3-3. The phase 
shifting frequency used in this simulation was fast enough 


to cause an increasing amplitude wander to be produced. 
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Figure 3-9 
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E. CONCLUSIONS 

Producing the roll wander with random noise signal 
ensured that the model was behaving properly. The model 
meecevea CO the introduction of noise in the same fashion 
as the real missile. 

The use of the harmonic generator to produce a sustained 
smooth roll wander proved that the phenomenon was in fact 
eeused by harmonic addition. In addition, the use of the 
harmonic generator gave valuable insight into the roll 
ieep's response to external stimulus. The fact that the 
output of the bistable element drives the control surfaces 
eeeouecn a high rate, producing such large missile roll rates, 
meorains why it is so difficult to view an imbalance in 


mie bistable output. 
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Vem CObDUTION 


fee)6LNTRODUCT ION 

Once the cause of the roll wander phenomenon had been 
isolated and verified through testing of the missile model 
merece NO1se Source and harmonic generator, the problem 
existed to find a way of preventing this wander. As dis- 
meased Carlier, two causes of wander due to a noise input 
exist. Measures to prevent large amplitude wander were not 
considered. Large amplitude wander would be encountered 
only under extremely abnormal conditions in which the noise 
peedal Lotally saturates the control loop. On the other hand, 
tow amplitude wander caused by harmonics present in the noise 
eaaqding with the fundamental dither frequency, was felt to 
Pemoeolenizticant problem. In this chapter, the term roll 
wander will be understood to deal with low amplitude wander. 

The most direct method of preventing roll wander would 
Semone way venUatton Of the dither harmonics in the roll con- 
trol loop. Piss was a deceptively difficult problem. There 
already existed in the control loop several high frequency 
filters to filter out those high frequency signals generated 
by the missile's elastic body bending modes. Notch filters 
vice a simple low pass filter had been utilized to prevent 
the shifting of the fundamental dither frequency due to the 
Phase woal1tb Amposedq by the addition of the filters. Inclu- 
sion of another set of filters notched at the main harmonic 


components of the dither fundamental could be used to combat 
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m@emtrequency addition and resulting roll wander. However, 
feepee the dithering action of the roll loop is dependent 

upon the phase crossover of the loop, which is determined by 
meow oop e€lements, the frequency of the dither is dependent 
on the roll loop parameters and varies from missile to 
missile. More importantly, the dither frequency varies 
eeenoly irom one set of flight conditions to another. Thus, 
pamemem filter set to filter out harmonics would of necessity 
have to have a notch wide enough for all frequency variations. 
iene turn, would cause an unacceptable decrease in the 
dither frequency. <A self-adaptive notch filter was thus 


proposed as a possible solution. 


B. NOTCH FILTER DESIGN 

A filter with conjugate imaginary zeroes located at the 
frequency of interest, the notch frequency, would provide 
total attenuation of the signal at that frequency. At fre- 
M@uenertes fiieher than the notch frequency, the filter would 
have a gain greater than unity. However, the filter must 
have. unity gain at all frequencies other than the notch 
ieedquemcy. To counter the gain of the filter produced by 
the conjugate imaginary zeroes, a double pole is required. 
For stability reasons, purely imaginary poles are not possi- 
ble. The poles must have an attendant term. A problem 
arises, however, when complex conjugate poles with non zero 
real parts are introduced to produce a realizable filter. 

There is a sign reversal at the notch frequency, which 


ipiiecomas lee wonase shitt about the notch. Since the zeroes 
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Seow pure amaginaries, having no real parts, there is no 
phase disturbance at other frequencies due to this shift. 
This is not the case with the complex poles. The value of 
SC in the denominator of the filter expression determines the 
amount of phase disturbance the filter produces. Frequency 
components at lower frequencies, specifically at the funda- 
mental firequency, receive a phase shift depending on the 
maemmrruce Of ££. Tne shift in phase produced by the filter 
ehanges the phase crossover point of the entire loop. The 
phase crossover sets the limit cycle oscillation (dither) 
ifmequency. The value of ¢ also determines the effective 
feeoviaieot Lne Tilvjer notch. An engineering trade-off is 
required, therefore, between the width of the notch and the 
change produced in the dither frequency. 

icy eos polesmst be added ay some higher ITrequency 
(B = 3000) to allow the response of the filter to roll off 
at high frequencies. 

imeem yodemene Was mace dn the design of this filter that 
the zeroes of the filter were to be pure imaginaries. This 
does not imply Uiatwconple. CONnjUgave Zeroes With non zero 
real parts would not also provide attenuation at the notch 
frequency. However, if the zeroes aiso had non zero real 
COMpEoOnents. tne .amouny Of attenuation about the notch fre- 
quency would be decreased. With the width of the notcn 
ciceoadveowcramincais taccor, and since imaginary zeroes are 
possible (in practice; zeroes with negligible real parts), 


the filter was designed as stated. Equation Cis is sche 
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PMoaieeoimeie tCransier funetion of the filter proposed. 


Ue Ss 2 (17) 
(Ss? 4 Ceo + Wo J ( S + 8 ) 


Figures 4-1 through 4-6 show the Bode magnitude and phase 
mme@cs  10r the filter's frequency response from 1 hz to 1.59 
iaze1Oor values of ¢ of 0.01, 0.1 and O.2 respectively. The 
magnitude plots show the notch formed at the center frequen- 
cy. (A frequency of 56 hz, twice the fundamental dither 
frequency was used.) From the phase plots can be seen the 
EmMeunt Of phase shift introduced by the filter. With ac 
of 0.2, a phase shift of -25° was experienced at 28 hz. 

Her © equal to 0.1, the phase shift was -17.7°; and with 2 
pewal to 0.01, the phase shift was -10.8°. 

The filter was tested outside of the missile model. It 
was found that adequate attenvation of a signal having a 
bandwidth of 0.2 hz centered at the filter notch frequency 
MoaoeeouloODtaimed With values of © less than 0.0/7. When the 
filter was tested in the missile model with ¢ set to 0.01, 
the frequency of the dither was reduced from 28.75 hz to 
26.3 hz; with tT set to 0.1, the frequency of the dither 
dropped to 24.5 hz, and with a ¢ of 0.2, the frequency 
duoppedmcOovapproximabely 22 hz. 

Meevas cdeciaea to uwse a value for ct of 0.1 for minimum 
disturbance in the loop with maximum attenuation of the 


harmonic frequencies. Two filters were inserted, in cascade, 
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Bode Magnitude Plot of the Notch Filter with ¢=0.01 
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FIGURE 4-2 


Bode Phase Plot of the Notch Filter 
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FIGURE 4-3 


Bode Magnitude Plot of the Notch Filter with 
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PEGURE 4-4 


Bode Pnase Plot of the Notch Filter with z=0.1 
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FIGURE 4-5 


DedewidGimuldesr lot Of the Notch Filter with-7=0.2 
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FIGURE 4-6 


Bogen enase Floc of the Notch Filter with r=0.2 
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Tie nNegehn frequencies set to filter out the second and 


third harmonic frequencies of the fundamental. 


C. SELF ADAPTIVE FEATURE 

Wme notch frequency of the filter is specified by the 
value of W 6 in Equation (16). The compromise concerning the 
feeue Of £, Which determines the width of the notch was 
based upon the assumption that a wold be within £071 hz 
Oi cme harmonic frequency to be filtered. However, as 
previously discussed, several factors, including the inclu- 
on OL ene filter itself into the loop, will cause the 
fundamental frequency to vary. When this happens, the 
bearmon1c Of the fundamental also varies, by a larger amount, 
meooors1Onave to the order of the harmonic in question. If 
the harmonic component moves in frequency outside the .2 hz 
feces tmovcen, the filter no longer has any effect. The filter 
must therefore be self-adaptive; it must have a variable 
notch frequency. 

The filter obtains its self-adaptive feature by first 
measuring the dither frequency, then setting the filter 
eOeenmetonus vOrproduce the correct notch. The filter pro- 
duces an initial transient response. Constantly changing 
tie wimalwper coelficients, which is effectively the same as 
Como tanitly eiserbime a diifereny filter, disturbs the control 
loop. This disturbance is seen as a random wander in the 
DeceaeoOn me lOopraaToulMmrnimize this effect, the filter coeffi- 


cients must be held constant or be required to change slowly 
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mpea smectbh fashion. This will only be the case if the 
Ceenicr remains fairly stable. 

PAPsnuciminebnc dlhumer frequency is produced in the 
Simulation because continuous integrations are being simu- 
lated by discrete numerical operations. Also, since the 
Seerut Of the bistable element is a pure Square wave, 
instantaneous frequency was most easily measured at this 
point. Any external influence which disturbs the relative 
pumzaesa Of the states of the bistable element would also 
mmeaguce a fluctuation in the instantaneously measured 
imeeguency. The measured value of instantaneous frequency 
peemnGu tnerefore remain constant. When the filter coeffi- 
cients were changed each integration interval, to correspond 
to this instantaneous frequency, a wander with a frequency 
Sree pproximavely 2.5 hz was observed. The filter was itself 
causing a roll wander phenomenon (Figure 4-7). 

iesesieeretiecy Was Saenificantly reduced by using an 
average value for the frequency. The frequency was measured 
each menace ion iivemucl mt Melmaverated OVer six cycles. 

As each new cycle was recorded, it was then averaged with 
the merous five cycles; thereby "fairing" out any frequency 
jitter. The result is seen in Figure 4-8. 

Mi MmODurOrcmunat in actual practice the construction of 
such a self-adaptive filter is not quite so Simple as the 
Simulation Of one. However, techniques do exist for the 


Pueowerreon Of suleh filters, either analog or digital. This 
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FIGURE 4-7 


fevso TIME with Self Adaptive Filter in Control Loop (no input) 
Meaptation using Instaneous Frequency 
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Note: Initial excursicn is due to the initial conditions 
used in the simulation. 

Note: Dither frequency is measured each integration interval. 
Filter coefficients are updated at the same rate. 
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irGUrReE, 4-3 


fe elle LemeocelieAdastive Filter in Control (no input) 
Adaptation uSing Averaged- Frequency 
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Note: Dither frequency is measured and averaged over 6 
cycles. Filter coefficients are updated using averaged 
value of frequency. 


TO 





@mesis concerned itself only with the filter's effect on 


er emeeconuro! oop, and not on filter construction. 


ieee ehoTING THE FILTER 

PoMcemru Wace che effect of the addition of harmonic 
frequencies with the fundamental dither frequency which 
produced the distorted waveform responsible for roll wander, 
both the harmonic generator and the filter were inserted 
into the loop just prior to the bistable element (Figure 
4-9). In the actual missile, it would be expected that such 
e@aiter would be placed at the output of the roll rate gyro; 
Since it is the roll rate gyro which inputs the actual 
Marmonic frequencies into the loop. 

Mbeineeiec tmmlwer imscvalled, the control loop model was 
Tested in the same fashion as described in Sections II and 
III. Figure 4-8 and Figures 4-10 through 4-14 show the 
Paoocnisenol tne Model fo no input, step input and a Tail 
anduced torque, respectively. ' These runs were all made with 
memial, tLdentacal- flignt conditions. The missile's response 
during high altitudes and speeds remained the same as before 


maclusi@m of the iilter. 


E. WANDER SUPPRESSION 

The runs just described showed that with the self- 
adaptive filter in the loop, the model still responded 
properly to input commands. The model was then tested again 
with the harmonic generator. Figure 4-15 and Figure 4-16 


clearly show the filter's effect. The roll wander is totally 
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FIGURE 4-10 


KF, vo Tite wrehn FLtter in Loop (no input) 
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FIGURE 4-11 


Be Se ibey  ch enn LOOpeands—-415° step input at t=3.0 sec 
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FIGURE 4-12 


K ,.F vs TIME with Filter in Loop and a -45° step input 
at £=3.0 sec. 
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PrEGURE 4-13 


‘vs TIME with 
ee t=—3,0 sec 


Filter in the Loop 
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ana a +3°- Step torque (60) 





FIGURE 4-14 


Mone ver MMi with Filter im the Loop and a +3° step torque (6. ) 
af f= 3.0 sec. 
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FIGURE 4-15 


Bombo e ge ethe Eitect or a Varying Phase Harmonic 
( Harmonic Commences at t=1.0 sec ) 
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BEGur: i116 


¢6/vs TINE with Filter in Loop and with a Varying Phase 
Harmonic ( Harmonic Commences at t=1.0 sec. ) 
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emcee ssoca by Une Novch filter. The transient disturbance 
seen in Figure 4-16 is caused by the initial "switching in" 
Of the harmonic signal at t equals one second. The initial 
peatece. Of the harmonic produces a change in the dither 
frequency, with an attendant change in the output signal of 
miewmarmOnic generator. Since the notch filter is using an 
averaged value of frequency in computing the filter coeffi- 
weembs, the sShiit in harmonic frequency is momentarily 
Ser roce Lhe range of the notch. Roll wander starts. As 
the averaged frequency approaches the actual loop frequency 
the notch shifts and the wander is suppressed by the filter. 
Figures 4-17 and 4-18 show the model's response to the 
feeoeelripuc. A noise signal will disturb the roll loop and 
eeoewune missile to roll. This, However, is not roll 
wemder., Comparing these two figures, it can be seen that 
the amplitude of the wander has been decreased by about 
50 per cent. The rougher appearance of the curve in Figure 
4-18 is due to the fact that the random disturbance of the 
noise is becoming more dominant as the harmonic addition is 


being suppressed. 


F. CONCLUSIONS 
The primary objectives of this thesis were: 
Mueecomsutuleuman accurate Gifital model of the missile 


etme © IO a LOOD. 


fe eOcdueco min tnesmodel the pnenomenon of roll wander 
with an explanation of the sources and mechanisms of 
the wander. 


3.20 produce @ means of preventing or at least signi- 
ficantly reducing the occurrence of this phenomenon. 
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FIGURE 4-17 


¢’ vs TIME Showing the Effect of Adding White Noise (constant 
Noise Commences at t=0.7 sec. | 
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FIGURE 4-18 


Mow liiPegitae es lter in Loop and White Noise Added 


Noise Commences at t 
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The first two objectives have been achieved in this 
Piero ns slow in sections 1, Il, and III, the model 
Poaieweorwim tne Same fashion as the real missile. The digital 
Samulation of the missile now exists at the Naval Postgraduate 
pemMool, and is available for further study. 

‘MeN: solution presented in this chapter as a cure to the 
roll wander problem works and will prevent this phenomenon. 
However, a self-adaptive filter has several weaknesses which 
Peo oecineculrveness. The inclusion of the filter into 
meewcontvrol roll has reduced the fundamental dither frequency 
Significantly. Such a reduction may cause interference 
er oorems Witn the missile guidance system. It was found that 
evoneing the filter coefficients, to change the notch fre- 
aie Glouurded the convrol loop, causing a fluctuation in 
the fundamental frequency and random wander in the position 
loop. To minimize this effect, an averaged frequency was 
(ieee becalce Of this, the filter did not respond to fre- 
Pienerschane@es as quickly; and the roll wander eaused by the 
noise input was ee totally suppressed. 

It Pomoc re veCduunav much room exists for future thesis 
work in the area of preventing the roll wander. The largest 
portion of the work involved in this thesis was the construc-— 
ploOnmar an a2ecurace model and the production of the roll 
wander. With the model now in existence, future thesis 
students can use it and concentrate on the development of 


vos ewe metidcemoany System to prevent roll wander. 


83 





VY. RECOMMENDATIONS FOR FUTURE STUDY 


Peer olscocann tie previous Chapter, continued study 
mowiececedq in the implementation of an adaptive filter for 
mre suppression of roll wander. The filter used in this 
thesis produced a decrease in dither frequency of almost 
4 hz. This alone may make the filter acceptable. 

mmevner area of study, apart from the study of the missile 
control system, also became evident during the course of this 
I[meois. That is the study and the use of digital simulation 
languages. As is often the case with any research, ie 
preparation consumed as much time as the investigation itself. 
iam referring to the production of the DSL plot package. 
DSL lay dormant at the Naval Postgraduate School because of 
Maomiieabiltpy LO Output to the CALCOMP plotter. Once an 
operational plot package was produced, the value and the 
power of DSL as a simulation tool rapidly became apparent. 

tees COmbanuous system Modeling Program (CSMP), another 
simulation Wapeaace at this school, also at present requires 
considerable programming to access the CALCOMP plotter. It 
wowla require vonly sient modification to DSLPLOT and minor 
procedural changes to the system control cards for CSMP to 
Make Come output to the plotter via DSLPLOT. 

It must be remembered that as with any simulation routine 
DSL is time consuming; and on a time-shared system such as 
the IBM 360/67 System used at the Postgraduate School, turn 


around time is, at best, marginal. The IBM 360/67 System 
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has interactive graphics terminals. If DSL were modified 
via appending FORTRAN programs to output and input via a 
Praphiiecs terminal, a pseudo real time system could be 
produced. The advantages of such a system would be great. 
the student could modify his system, and automatically see 
the result. If the program misbehaved, the student could 
menace the run immediately, rather than waiting the full 
meceii led time of the run and the time required to exit the 
eyoavem prior to being able to study the program and take 
corrective action. sSuch a system could easily allow the 
Mocteue accomplish in a two or three hour session what might 
require days or weeks to accomplish off line. The modifica- 
feemeol DSL to produce such a system would be extensive but 
Rewarding. 

The other modeling languages on file at the Postgraduate 
pemool (CSMP, ECAP, LISA, etc.) all would benefit immensely 
eaeowmilar modifications. The saving in systems time by 
removing such jobs from the normal job queues and rescheduling 
Maem. Gurins hours of minimum load would also be significant. 

The emphasis in control engineering has shifted dramati- 
cally to Bae Miz aipmeme Ot tne digital computer, both for 
Prootem solving and system modeling. However, the structure 
and implementation of the simulation programs available is 
Heumexoliciuly taucht., A course bULItT around system modeling 


using such programs would be highly beneficial. 
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APPENDIX A 


i INTRODUCTION 

Digital Simulation Language (DSL) is an IBM System/360 
MeminaNn LV program for the digital simulation of continucus 
System dynamics. Its nonprocedural problem-oriented input 
language accepts problems expressed either at the analog 
block diagram level or as systems of ordinary differential 
Povlavions. osoince the input language is nonprocedural, proper 
maravemeny Sequence is established by an internal sort based 
Semeecorrecy antormation flow. DSL includes FORTRAN as a 
Subset, thereby extending its power to handle nonlinear and 
time-variant problems of considerable complexity. <A trans- 
fewer converts the input statements into a FORTRAN subprcgran, 
Yoel 1S compiled and executed to accomplish the simulation. 

A significant feature of DSL is centralized integration, 
which means that all integrator outputs are computed simul- 
wameously av the end of each calculation interval. The user 
may specify one of several integration methods, including 
fifth order Milne Predictor and Corrector and second and 
mourtn orcger Runga Kutta. 

Dsl as Originally presented at the Naval Postgraduate 
SSNeoGimoulupuus, vo the high speed printer only, either as 
Dicmvead GUTDUG Or printer plots. Due to software incompat— 
tPivicc we Dolmas ailcrally installed would not output to 
the CALCOMP XY plotter. This thesis, as with many simulation 


Studies, required the CALCOMP as an output. For this reason, 
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and since DSL had so many advantages, extensive work was 

Some prior LO Constructing the missile simulation on a FORTRAN 
plotting program to output the DSL program via the CALCOMP 
meovover., Inis program, DSLPLOT and its associated FORTRAN 


peorram PLOTHALL, will be discussed later. 


B. PROGRAM FLOW 

emcee Start Of the program, all program variables are 
ferme oeeme initial conditions given. TIME, the independent 
variable, is set to zero, and a run counter KSIM is incre-~ 
mented if it is the beginning of a new parameter study. 
mreally KoliM 35 sev equal to 1. At this point, the region 
Veamewprorran Called the INITIAL Region is entered. In the 
INITIAL Region, arithmetic operations and logical statements 
to be executed only at the beginning of a particular run are 
Comeuuved. An integer flag KEEP, is set equal to zero, and 
the DERIVATIVE Region is called. 

The DERIVATIVE Region represents that portion of the 
muha on wheen involves integration and the calculation of 
the derivatives. The basic interval in the independent 
variable for each pass through this region is the calculation 
tomenval DELI. It 1s for this reason that the DERIVATIVE 
Region should include all time (independent variable) 
Gependent calculations which must be performed at the same 
frequency required by the selected integration algorithm to 


assure accurate numerical integration. 
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iiemicmreocedmral Navure Of Lhe input language is the 
GeetmipeO: bhe DERIVATIVE, Region. The DERIVATIVE Region is 
Moiomeoney  —TreeTon in whieh an intermal sort routine is utilized 
tO establish the correct statement sequence. Only structure 
Stacvements are sorted and sequenced. A statement is properly 
eeeuenced if all its inputs are available either as input 
meepemevers OF initial conditions or as previously computed 
Values in the current calculation interval. Outputs of INTGRL 
and MEMORY blocks are known because their initial values are 
Biven. The internal sequencing procedure begins first with 
icuerralOr OuUupuLSs, then with outputs of MEMORY blocks and 
tonoally with Chose remaining output variables in the 
Peatveli ve Region that have not been ineluded in the sort 
sequences of integrators and MEMORY blocks. It is important 
Wemirecelize that the sorving algorithm is unable to sequence 
the statements properly unless it can assume a known output 
variable as a starting point of a sort sequence. AS previously 
stated, this output is provided by either an INTGRL or MEMORY 
block. A loop which does not contain at least one such block 
iomcowmuea ran EMiPLiCiT or ALGEBRAIC Loop and cannot be 
Sequenced..) Upon encounverings Such a loop, the sorting algo- 
rithm is aborted, and all subsequent statements are executed 
Peon Ge Ir. 

The use of the sort algorithm places a severe restriction 
on the type of statements which can be used in the DERIVATIVE 
Region. No FORTRAN logic or branch statements can be used 


because they cannot be sorted properly. To allow the 
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versatility which accompanies these types of statements, 
emeenwecvery incvegrarion interval, the program flow exits the 
DERIVATIVE Region and enters the DYNAMIC Region. The DYNAMIC 
Peron provides an area in which FORTRAN logic can be used. 
Hestins Of variables, adjustment of coefficients and nonlinear 
@eetatton all can be accomplished. It should be remembered, 
memever, that since for a given value of the variable TIME, 
the DERIVATIVE Region is executed before the statements in 
the DYNAMIC Region a time lag equal to DELT exists between 
Presacvuions taken in the DYNAMIC Region and the resulting 
eeieey an the DERIVATIVE Region. 

Pececona Sysvem Variable, DELS, dictates the frequency 
Seevoien the computation interval is interrupted and program 
meowels diverted to the SAMPLE Region. The makeup of the 
SAMPLE Region is the same as that of the DYNAMIC Region. 
feewinaim Lunetion iS to test variables to determine points 
for data collection and to order the program flow to allow 
Simultaneous collection of data of several variables. 
DSLPLOT, which wae iiovembedeavstme Naval POsiteraduatre Ocnool, 
is accessed via the SAMPLE Region. 

Finally, at the termination of the simulation either by a 
logic statement or when the variable TIME is set equal to the 
system variable FINTIM (a variable set by the user), the 
program makes one last pass through the SAMPLE Region and 
enters the TERMINAL Region. In the TERMINAL Region, logic 
cam oeemlised TOeMOdiiy Orogram variables and constants and a 


new run initiated. If DSLPLOT has been utilized, the plotting 
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Sevomeer to Ccloseq Out. If no action is taken to initiate 
pemewerum, the DSL job will be terminated. If a new run 
mes Deen Called for program flow returns to the INITIAL 
Region. 

When the DSL job is completed, the DSL program is 
dynamically linked via the IBM 360 linkage editor to the 


erogzram DSLPLOT. 


¢ 


fee DoLPLOT 

The program DSLPLOT was produced as a forerunner to this 
thesis to allow the DSL output variables to be displayed via 
meee CALCOMP AY Plotter. DSLPLOT accomplishes this by writing 
tne variables to be plotted along with counters indicating 
Wome rou nMMumMoer., curve number, and individual point number on 
eumexvernal data set via the subroutine DRWG each sample 
interval (DELS). Each group of data points indicating a 
Single DSL run is delineated by a set of zeroes produced by 
the subroutine ENDRW called in the TERMINAL Region. 

PLOTHALL, the. appending FORTRAN program used by DSLPLOT, 
ereecesses Uhis data set Upon a successful completion of the 
DSL job. A three-way sort is performed on the data set, and 
mie pOlnts to be plotted are written in order on a second 
data set. The points are then retrieved one curve at a time 
from the second data set and passed to the plotting subroutine. 
This subroutine utilizes the IBM Plotting Package and produces 


the required CALCOMP plots. 
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since it was intended that DSLPLOT was to be a permanent 
eoemo1On tO the DSL program at the Naval Postgraduate School, 
meewes Written With as much flexibility as possible. Listed 
i@eeleote A-l are the various modes under which DSLPLOT can 
operate. The different modes allow the user to choose the 
size of the output plot and to either specify the scaling 
used on either axis or to allow the program to auto scale one 


or both axes. 


Peeewoe OF THE PLOT PACKAGE 
MoeprevTouwsly discussed, points are collected for plotting 
via the subroutine DRWG. A call to this subroutine is 


Haeeraved by Che FORTRAN statement: 


CALL DRWG(1,1, NUM,X,Y) 


Where X and Y¥ are the variables to be plotted. The xX 
Variable is plotted as the abscissa coordinate and Y the 
ordinate value. NUM is the integer counter denoting the data 
point number. The two integer constants appearing as the 
farst two arguments represent the plot number and curve 
number for that plot number. These arguments may be provided 
as integer variables. This allows overlaying a set of curves 
from different runs on the same plot. It must be remembered 
Dimer aoremrs No Comvention in DSL, as there is in standard 
HORTRAN, for providing integer variables. Those variables 
such as NUM, or the variable NPLOT used in ENDRW must be 


declared integers by a DSL INTGER statement. 
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TABLE A-1 


values auto scales will be produced. 
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Mode No. User Supply (in order listed) 
0 Auto (no inputs required) 
1 Pinte > eee Int Dy 
2 Amin Dx 
S Ymin Dy 
4 Lenx Leny 
5 Xmin Dx Ymin Dy lLenx Leny 
6 Minin Dx Lenx Leny 
7 Ymin Dy Lenx Leny 
8 Lenx 
9 Xmin Dx Ymin Dy Lenx 
10 Amin Dx Lenx 
dla Ymin Dy Lenx 
le Leny 
lag) LSyabiel Bes Leny 
14 Xmin Dx Ymin Dy Leny 
3) Ymin Dy Leny 

Note: If mode 0 is specified a 9" by 9" graph with all 





Upon completion of data collection for a given set of 


merous, vene data set is closed out by the FORTRAN statement: 
CALL ENDRW(NPLOT) 


ftietemirhOl is an inveger variable which is initially set to 
the number of intended data sets. Each call to ENDRW 
decrements the variable by one. When NPLOT is less than or 
eaelepe Zero, the DSL job 1s terminated. After a call to 
ENDRW, if NPLOT is still greater than zero, the program 
Simply exits the subroutine, and the DSL program continues. 

Bach data set must be closed out by a call to ENDRW. 

If the plotting package has been used in the job, ENDRW must 
aiso be used as the means of terminating the DSL job. 
PEOTHALIL, link edited to the DSL job, terminates upon encoun- 
veri aevalue= on NPLOT equal to one, therefore if the main 
Peon 1S) cverminactved by means other than ENDRW, an abnormal 
termination of the plot package will result. 

Since it is standard practice to initiate a new set of 
plots each DSL run, the plot, curve and data point indicators 
refer to a given run. After a call to ENDRW, if a new run is 
initiated either as a parameter study or by a call to RERUN, 
CONTIN, etc., encountering the same DRWG statement again 
will generate a new plot. If curves are to be overlaid on 
one plot; at the end of a given run the data set must be held 
open by not calling ENDRW and the curve indicator, specified 


by a variable name, is updated. When the plot is complete, 


e' 





moe wenwoct then be called to close the data set. It is 
mipervant tO nove that in any case, NUM, the data point 
Seiicer must be initialized to 0 at the beginning of each 
Peaeanceiieremenved prior tO its respective group of DRWG 
BeaeemMents . 

It is important-also to remember that ENDRW merely closes out 
Macmcurrent plot data set and decrements its calling argu- 
fete WOOL. ENDRW does not initiate a new DSL run. This 
fueteoc done by the user with a call to RERUN, or CONTIN or 
Dy a Darameter study or by any other method given in the 
DSL USER'S MANUAL. Figure A-1l is a sample outline of a DSL 
Megeran, Snowing the use of the subroutines DRWG and ENDRW 
BOeproduce plot output. In Figure A-l1, two data sets are to 
bemmormed. Hach data set will contain two plots. The first 
Prous iINeeacn Gata sev will contain two curves. 

Another important point is that the subroutines DRWG and 
ENDRW are completely independent of the DSL statements 
PREPAR, GRAPH and PRTPLOT. The original IBM system output 
its data to the AY plotter via the statements GRAPH and 
PREPAR. At the Naval Postgraduate School, however, PREPAR, 
GRAPH and PRTPLOT (which at this installation must accompany 
GRAPH) are only used to produce the high-speed printer plots. 

The plot package will handle up to ten plots per run and 
ven curves per plot. Any given curve can have up to 2000 
Eometso mI anorescnad c©O0O0 perms are collected for a curve, 


the DSL run will be terminated, and the TERMINAL Region entered. 
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INTGER NUM,NPLOT 

CONST NPLOT=2 

CONTRL FINTIM=10.0,DELT=0.1,DELS5=0.1 
INITIAL 


NUM=0 
DERIVATIVE 


SAMPLE 
NUM=NUM+1 
CALL DRWG(1,1,NUM, TIME, Xx) 
CALL DRWG(1,1,NUM,TIME,Y) 
CALL DRWG(2,1,NUM,X,Y) 


TERMINAL 
CALL ENDRW(NPLOT) 
(parameter or constant changes 
using either FORTRAN assignment 
statements or a DSL END statement 


' followed by the appropriate DSL 
non-executable statements) 


CALL RERUN (when FORTRAN assignments 
have been used) 


END 
STO 


Figure A-1 


oe) 





leeienmcceosmuacomouletron Of the DSL job, the plot 
eaorran PEOTHALL is executed. Each plot set up in the DSL 
HOb Must be accompanied by two title cards. Columns 1 
through 48 on both cards are reserved for titles. Columns 
75 through 80 of the second title card constitutes the mode 
field (Figure A-2). 

The mode numbers are listed in Table A-l. If a printed 
SNtouL 1S Specified by supplying a printed output interval, 
then a listing of the points actually plotted will be provided 
fooled watetnemlnverval sSpeciiied. The printed output 
interval can have a range from every point plotted to one of 
every 999 points plotted (assuming 999 or more points exist). 
Pacoolumme 7) 1s lelt blank, no grid of scaled axes will be 
drawn; just a left-hand and bottom scale will be provided. 
Cimawieroeolaced an this column, a labeled grid at 3 inch 
iitcervals Will be drawn. This grid is time consuming and 
PneuGmeoe ONLtTeoed 1f Met essential. 

If a mode other than mode 0 is specified a third title 
card must be provided hOmethie Dartcular ploy in question. 
The appropriate data must be entered via this title card 
as eho aii Figure A-3. 

Figure A-4 provides a listing of the JCL (job control 


language) cards necessary to implement a DSL job. 
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Columns 1 thru 48 on the 
ilmeespesion cle le, cards per 
POM muemer bln le 11e1d . 
The user's name or box 
number shoulda appear 
somewhere within this field 
Omeeltner card. 








Mode Field: 
Column 75: ~erid daesired 
Column 76-78: printed output 


iinwetwe 
Column 79-80: mode number 


Figure A-2 
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M4 
Scaling Fields Graph Dimension Fields 


Note: Columns 1 thru 40 are reserved for Xmin, Dx, 
Ymin and Dy (fields left justified as needed). Values 
ower tiven eT ormat Or & format. if E format is 
MecOmmecaweiiioceDe TIphne JUSCiiled in the appropriate 
field. 


Note: xXmin is the value appearing on the extreme left 
SMOmOoumurennorlzZontal scale @nd Ymin is the value 

reo omanewena Gae DobuCm Of Vine vertical scale. Dx and 
tyeorewe mee tesuccurvye Scale increments. Lenx and Leny 
CRemememermensleons. of the plct in inches in the £ and 

facireetrons, where: 


Q< Lenx < 9.0 
QO < Leny < 20.0 


Seviousty hx and Dy must be compatible with Lenx and 
Leny. When no lengths are specified a 9.0" by 9.0" 
Plotewtilt be drawn. 


Note: Columns 41 thru 45 and 51 thru 55 are reserved 
for Lenx and Leny (fields left justified as needed). 
VerUesmmuss Oe anpuc in F format. 


Figure A-3 
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V7 (standard green job card) 
// EXEC DSL 
7705 ae NeUT DD: * 


(dsl deck goes here) 


V7 eVOleeototie DD * (if plots are to be drawn) 


(2 or 3 title cards per plot requested) 


/* 


jk (Note: 2 /*'s must be provided) 


Figure A-4 
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